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Abstract 

We show how the entropy operators for two subsystems may be calculated. In the case of the 
atom-field interaction we obtain the associated Wigner function for the entropy operator for the 
quantized field. 

PACS numbers: 02.10.Ud, 42.50.Dv, 03.65.Ud, 42.50.Lc 



1 



I. INTRODUCTION 



It is well known that when two subsystems A and B interact the entropy of the total 
system and the subsystem entropies obey a triangle inequality [1] 

\Sa — S B \ < Sab < Sa + Sb- (1) 

This result implies that if both subsystems are initially in pure states, the total entropy is 
zero and both subsystems entropies will be equal after they interact. Here we would like to 
arrive to the result that, if initially the two subsystems are in pure states any function the 
density matrix of subsystem A is equal to the function of the density matrix of subsystem 
B, and in particular obtain that both subsystems entropies are equal, without using the 
Araki-Lieb theorem. We would also like to find the entropy operator 

Sa(b) — — InpA(B), (2) 

for any of the subsystems. In particular we will consider later a two-level atom-field interac- 
tion, but the results may be generalized to other kind of subsystems, for instance atom-atom 
interaction, atom-many atoms interaction, iV-level atom-field interaction, etc. 

The density matrix for a two-level system interacting with another subsystem B is given 

by 



I \c){c\ 


\c)(s\ \ 


\\s)(c\ 


\s)(s\ 



P= (3) 



where |c) (\s)) is the unnormalized wave function of the second system corresponding to the 
excited (ground) state of the two level system. From the total density matrix we may obtain 
the subsystem density matrices as 

^ =| W(*)\ (Vp 12| _ (4) 

x (C|S) {S\S} J yp 2 l p22 

and 

Pb = \c){c\ + \s)(s\. (5) 

From the density matrix for subsystem B one can not see a clear way to calculate the entropy 
operator as powers of ps get complicated to be obtained. To make the calculation easier we 
state that if two subsystems are initially in pure states, after interaction, the trace of any 
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function of one of the subsystems density matrix is equal to the trace of the function of the 
other subsystem's density matrix. 

We can prove this, by using the following relation valid for two interacting systems that 
before interaction were in pure states 



PB + =Tr A {p(t)p n A (t)}, 



with this relation it is easy to show that 



Tr B {pT(t)} = Tr A { P n A +L (t)} 



(6) 



(7) 



In particular with the expression (6) we demonstrate that Sb — S A . Of course it is also true 
that 

p n A +1 =Tr B {p(t)p B (t)}. (8) 
II. PROVING pl +1 = Tr A { P (t)p n A {t)} BY INDUCTION 



We can prove relation (6) by induction. For this we need to find p A , to this end we write 



Pa 



9 Pl2 



+ R, 



(9) 



\p21 -J 

where S — pu — P22 and 1 is the 2x2 unit density matrix. We can then find simply 



Pa 



n 



-BP 



(10) 



=0 \ m 



We split the above sum into two sums, one with odd powers of R and one with even 
powers. We also use that 



^2m ^2m 



R 2m+1 = -6 



2m+l 



/ 2 \l/2 

with e = I V + I 

P12I I • Therefore we can write 



« = 5 



, R 

1+ 2? 



1 

2" £ 



In terms of p A the above equation is written as 



(11) 



(12) 



p\ = G(n)p A - \\p A \\G(n- 1)1 



(13) 
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where 

G(») = h 

with the determinant ||pa(£)|| — \ — ^ ■ Note that we have written p\ in terms of p A and 
the unity matrix. We could have done it via Cayley-Hamilton theorem [2], which states that 
any (square) matrix obeys its characteristic equation, i.e. a 2 x 2 matrix may be written as 
we did, as the square of the matrix (and any other power of it) may be related to the matrix 
and the unity matrix. 

Then relation (6) my be proved by induction: for n — 1, we have 

p 2 B = Tr A {p(t)p A (t)}. (15) 

We now assume it to be correct for n = k + 1 

p k B +1 = Tr A {p{t)p k A {t)} (16) 

and prove (6) for n — k + 1. By using (13) we can write 

p k B +1 = p%G(k) - p B \\p A (t)\\G(k - 1) (17) 

Note that any power of p B may be written in terms of p B and p\ (k > 2). Multiplying the 
above equation by pB we obtain 

P B +2 =plG(k)-pl\\p A (t)\\G(k-l) (18) 

We obtain p 3 B from (6) as 

P% = Tr A {p(t)pl(t)} = Tr A {p(t)lPA(t) + (e 2 - \)}} (19) 

where for the second equality we have used the Cayley-Hamilton theorem for the atomic 
density matrix: p A (t) = p A (t) — Inserting (19) in (18) and after some algebra we 

find 




PT = PlG{k + 1) - p B \\p A (t)\\G(k) (20) 



or 



Ps + =Tr A {p(t)PA +1 (t)} (21) 
that ends the prove of relation (6) by induction. 



4 



III. ATOMIC ENTROPY OPERATOR 



With the tools we have developed up to here we can study the two-level atom-field 
interaction [3] or equivalently the ion-laser interaction [4, 5] and construct atom and field 
entropy operators. In the off-resonant atom-field interaction, i.e. the dispersive interaction, 
the (unnormalized) states |c) and \s) read [6] 

\c) = ^=\(3e-^), \8) = ±\f3f*) (22) 

where \/3) = e - ^ 2 / 2 YlkLo 7^1^) * s ^ ne hiitial coherent state for the field, x is the interaction 
constant. The atom was considered initially in a superposition of excited and ground \ipA) = 

^(|e> + b». 

We write the entropy operator as 

S A = lnp^ 1 = ln(l - p A ) - In \\pa\\, (23) 

such that the expectation value of Sa is the entropy. In the above expression we have 
used that p^ 4 1 (t) = £a/||pa(£)|| with the purity operator £a = 1 — p A (t). We can develop 
ln(l — x) = — Yl^Li ^T an< i use (13) to find 



S A = F 1 p A + F 2 1 (24) 



with 



F 2 = — 
2 



2e V 1 + 2e 
1 



\n\\p A (t)\\ + -\n 



2e V 1 + 2e 



(25) 



Note that Sa is linear in the atomic density matrix as expected from Cayley-Hamilton's 
theorem (for 2x2 matrices). From (24) we can calculate the atomic (field) entropy and the 
atomic (field) entropy fluctuations 

(S A ) = F 2 + F 1 (l-2\\p A (t)\\), (26) 

and 

(AS a) = y/(S\) - (S A ) 2 = In IM*)I| 1/2 - (27) 
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IV. FIELD ENTROPY OPERATOR 



We use the expression for p B in terms of p A to write the field entropy operator in terms 
of the atomic density operator 

S B = Tr A {p{t)S A {t)p A 1 {t)} (28) 

Note that we can write the atomic entropy operator in terms of the operator used to define 
concurrence [7] because p~ A \t) = p^^P^(0^ = ]0$\\ i- e - 

Sb = -—±-^Tr A {p(t)S A p A (t)} (29) 
\\Pa( 



By inserting (24) into (28) we obtain 

S B = Tr A {p(t)(F 1 + F 2 p A 1 (t))} (30) 

using the expression of the inverse of the atomic density operator in terms of the purity 
operator,the entropy may be written as 

S B = Tr A {p(t)(Fi + jT^^ 1 " M*)])} (31) 
In terms of the field density matrix the field entropy operator is 

Sb = (f, + — ^) p B (t) - jr^wrMt) (32) 

V \\PA{t)\\J \\PA{t)\\ 

from (5) we can write p B (t) as 

p 2 B = \c)(c\(c\c) + \s)(s\(s\s) + \c)(s\(c\s) + \s)(c\(s\c). (33) 

From (22) we have that (c|c) = (s\s) = 1/2 and that (c\s) = (s\c)* = exp(-|/3| 2 [l -e 2ixt ])/2. 
The entropy operator for the field is written as 

Sb = U + ^twu) (\c)(c\ + \s)(s\) - TT^wn((s\c)\s)(c\ + (c\s)\c)(s\) (34) 

We have chosen the dispersive interaction so that all the terms forming the field density 
matrix and its square are coherent states, such that we can calculate the Wigner function 
associated to (32) in an easy way. We do it by means of the formula[8] 

oo 

W s (a) = ^2(-l) n {a,n\S B \a,n) (35) 

n=0 
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where \a,n) are the so-called displaced number states [9]. The explicit expression for the 
above equation is 

W s (a) = e -^-2|aP + 4/^cosx* (p, + 2|| ^ t) J cosh(4/^sin X t) 



e ~2^~2\a\*+4fia x cos X t F2 cos(2/3[/3 sin 2 X t - 2a x sin x*])- (36) 

2 \\PA{t)\\ 



V. CONCLUSIONS 



We have shown that the entropy operator for the subsystems may be obtained and in 
the case one of those subsystems is the quantized field. The conection with quasiprobability 
distribution functions [10, 11], in particular the association to the Wigner function for the 
entropy has been obtained. 
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